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We reconstruct B — > DK ± decays in a data sample collected by the CDF II detector at the 
Tevatron collider corresponding to 1 fb~ of integrated luminosity. We select decay modes where 
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and measure the direct CP asymmetry Acp+ = 0.39 ± 0.17(stat) ± 0.04(syst), and the double ratio 
of CP-even to flavor eigenstate branching fractions Rcp+ = 1.30 ± 0.24(stat) ± 0.12(syst). These 
measurements will improve the determination of the Cabibbo-Kobayashi-Maskawa angle 7. They 
are performed here for the first time using data from hadron collisions. 
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The measurement of CP asymmetries and branch- 
ing ratios of B~ — > DK~ [l| decay modes allows a 
theoretically-clean extraction of the CKM angle 7 = 
arg(— V u dV* b /V c dV* b ), a fundamental parameter of the 
standard model [2(. In these decays the interference be- 
tween the tree amplitudes of the b — > cus and b — > ucs 
processes leads to observables that depend on their rela- 
tive weak phase (7), their relative strong phase (5b), and 



the magnitude ratio tb 



A(b- 



A{b- 



These quantities can 



all be extracted from data by combining several experi- 
mental observables. This can be achieved in several ways, 
from a variety of D decay channels 0-01 • 

An accurate knowledge of the value of 7 is instrumental 
in establishing the possible presence of additional non- 
standard model CP-violating phases in higher-order dia- 
grams @, 0] ■ Its current determination is based on a com- 
bination of several B — > DK measurements performed in 
e + e~ collisions at the T(4S) resonance [cHlPf and its un- 
certainty is between 12 and 30 deg, depending on the 
method [ll[. This uncertainty is almost completely de- 
termined by the limited size of the data samples avail- 
able, with theoretical uncertainties playing a negligible 
role (~ 1%). The large production of B mesons available 
at hadron colliders could offer a unique opportunity to 
improve the current experimental determination of the 
angle 7. However, the feasibility of this kind of mea- 
surement in the larger background conditions of hadronic 
collisions has never been demonstrated. 

In this paper we describe the first measurement of 
the branching fraction ratios and CP asymmetries of 
B~ — > DK ~ modes performed in hadron collisions, based 
on an integrated luminosity of 1 fb _1 of pp collisions at 
y/s = 1.96 TeV collected by the upgraded Collider Detec- 
tor (CDF II) at the Fermilab Tevatron. We reconstruct 
events where the D meson decays to the flavor-specific 
mode K~ir + (D®), or to one of the CP-even modes 

K~K+ and 7^+ [D C p+ = (D°+S°)/V2|. From these 
modes, the following observables can be defined: 



.4 



CP+ 



B(B~ 



D 



CP+ 



K-)-B{B- 



D 



CP+ 



K+) 



B(B~ -> Dcp+K-) + B(B+ -> D C p+K+) ' 



(1) 
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CP+ 



B(B- 



D 



CP+ 



K- 



B(B^ 



D C p+K+) 



B(B- -> D° f K-)+B(B+ D f K+) 



(2) 

With the assumption of no CP violation in D° decays, 



and neglecting D°— D° mixing 12 1. these quantities are 
related to the CKM angle 7 by the equations Q 



Rcp+ = 1 



2r cos 8b cos 7, 



A C p+ = 2r B sin i5b shi7/P C p + . 



(3) 
(4) 



For our measurements we adopt the usual approx- 



imation Rcp+ ' 
r ■ \V us V cd /V ud V ct 

B{B~ 



R 



j-, which is valid up to a term 
0.01 fil, where 



R = 



D° f K-)+B(B^ 



D° f K+) 



R4 



B{B~ 



Dcp+K-) + B{B+ -> D C p+K+) 



B(B~ -> Dop+tt-) + B(B+ -> Dcp+tt+) 



(5) 
(6) 



The CDF II detector is a multipurpose magnetic spec- 
trometer surrounded by calorimeters and muon detec- 
tors. The components relevant for this analysis are briefly 
described here. A more detailed description can be found 
elsewhere [l4|. Silicon microstrip detectors (SVX II and 
ISL) [H and a cylindrical drift chamber (COT) (HI im- 
mersed in a 1.4 T solenoidal magnetic field allow recon- 
struction of charged particles in the pseudorapidity range 
I 77 |< 1.0 [13 • The SVX II detector consists of mi- 
crostrip sensors arranged in five concentric layers with 
radii between 2.5 and 10.6 cm, divided into three con- 
tiguous sections along the beam direction z, for a total 
length of 90 cm. The two additional silicon layers of 
the ISL help to link tracks in the COT to hits in the 
SVX II. The COT has 96 measurement layers between 
40 and 137 cm in radius, organized into alternating ax- 
ial and ±2° stereo superlayers, and provides a resolu- 
tion on the transverse momentum of charged particles 
a p T IVT — 0.15%pr/(GeV/c). The specific energy loss 
by ionization (dE / dx) of charged particles in the COT 
can be measured from the collected charge, which is en- 
coded in the output pulse width of each sense wire. 

Candidate events for this analysis are selected by a 
three-level trigger system. At level 1, charged par- 
ticles are reconstructed in the COT axial superlayers 
by a hardware processor, the extremely fast tracker 
(XFT) [lj|. Two oppositely charged particles are re- 
quired, with transverse momenta pp > 2 GeV/c and 
scalar sum pp\ +pp2 > 5.5 GeV/c. At level 2, the silicon 
vertex trigger (SVT) [lj| associates SVX II r — <f> position 
measurements with XFT tracks. This provides a precise 
measurement of the track impact parameter, do, which is 
defined as the distance of closest approach to the beam 
line. The resolution of the impact parameter measure- 
ment is 50 /Ltm for particles with pp of about 2 GeV/c, 
including a w 30 /im contribution due to the transverse 
beam size, and improves for higher transverse momenta. 
We select B hadron candidates by requiring two SVT 
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tracks with 120 < e?o < 1000 fxm. To reduce background 
from light-quark jet pairs, the two trigger tracks are re- 
quired to have an opening angle in the transverse plane 
2° < Acf> < 90°, and to satisfy the requirement L xy > 200 
/jm, where L xy is defined as the distance in the transverse 
plane from the beam line to the two-track vertex, pro- 
jected onto the two-track momentum vector. The level 1 
and 2 trigger requirements are then confirmed at trigger 
level 3, where the event is fully reconstructed. 

Reconstruction of B~ hadrons begins by looking for a 
track pair that is compatible with a D° decay. The in- 
variant mass (Md) of the pair is required to be close to 
the nominal D° mass (1.8 < M D < 1.92 GeV/c 2 ). This is 
checked separately for each of the four possible mass as- 
signments to the two outgoing particles: K +/ n~ , A'~7r + , 
K + K~ and -k + it~ . The D° candidate is combined with 
a negative charged track in the event with pt > 0.4 GeV 
to form B~ candidates. A kinematic fit of the decay is 
performed by constraining the two tracks forming the D 
candidate to a common vertex and to the nominal D° 
mass, the D candidate and the remaining track to a sep- 
arate vertex, and the reconstructed momentum of the 
B~ candidate to point back to the luminous region in 
the transverse plane. 

To complete the selection, further requirements are ap- 
plied on additional observables: the impact parameter 
(<1b) of the reconstructed B candidate relative to the 
beamline; the isolation of the B candidate {Ib) the 
goodness of fit of the decay vertex (xb)', the transverse 
distance of the D, both relative to the beam \L xy {D)] 
and to the B vertex [L xy s{D)], and the significance of 
the B hadron decay length [L xy (B) j 'a i (b)]. We chose 
the requirement L xv b(D) > 100 [im to reduce contami- 
nation from (nonresonant) three-body decays of the type 
B + — > h + h~h + (from here on, we will use h to indicate 
either K or w), in which all tracks come from a common 
decay vertex. In addition, we reject all candidates com- 
prising a pair of tracks with an invariant mass compatible 
with a J/ip — > decay within 2a. The threshold val- 

ues for all other requirements, whose purpose is to reduce 
combinatorial background, were determined by an unbi- 
ased optimization procedure aimed at achieving the best 
resolution on Acp+- This resolution was parametrized as 
a function of the expected signal yield S and background 
level B, by performing repeated fits on samples of sim- 
ulated data extracted from the same multidimensional 
distribution used as likelihood function in the fit [Eq. 
[7]. For each choice of thresholds, the signal S was deter- 
mined by rescaling the number of observed B~ — > Z3y7r~, 
and the background B was determined from the upper 
mass sidebands of each data sample (5.4 < Mb < 5.8 
GeV/c 2 ). Based on this optimization procedure, we 
adopted the following set of requirements: Ib > 0.65, 
X 2 b < 13, d B < 70 /an, L xy {B)/a Ln:y{ B) > 12, and 
L xy (D) > 400 urn. 

For every B~ — > Dhr candidate, a nominal invariant 
mass is evaluated by assigning the charged pion mass to 



the particle h~ coming from the B decay. The distri- 
butions obtained for the three modes of interest (D — > 
Kit, KK or irir) are reported in Fig. [TJ A clear B~ — > 
Dtt~ signal is seen in each. Events from B~ — > DK~ 
decays are expected to form much smaller and wider 
peaks in these plots, located about 50 MeV/c 2 below the 
B~ — > Dtt~ peaks, and as such cannot be resolved. The 
dominant residual backgrounds are random track combi- 
nations that meet the selection requirements (combina- 
torial background), misreconstructed physics background 
such as B~ -> D*°tt~ decay, and, in the D° — > KK fi- 
nal state, the nonresonant B~ — > K + K~K~ decay, as 
determined by a study performed on CDF simulation. 

We used an unbinncd likelihood fit, exploiting kine- 
matic and particle identification information from the 
measurement of dE/dx in a similar way to [2lj . to sepa- 
rate statistically the B~ — > DK~ contributions from the 
B~ — > Dtt~ signals and from the combinatorial back- 
ground. To make best use of the available information, 
wc fit the three modes simultaneously using a single like- 
lihood function, to take advantage of the presence of pa- 
rameters common to the three modes. 

The likelihood function is 



PID 



PID 



(7) 



where c labels combinatorial background quantities, b 
is the combinatorial background fraction, and £ km and 
£ PID are defined below. The index j runs over the 
modes B~ — > DK~, B~ — > Dtt~ , nonresonant B~ — > 
K+K-R- and B~ n + ir~K-, and B~ ->■ D*°ir- 
(wherc a soft 7 or ir° from the D*° is undetected) and 
fj are the fractions to be determined by the fit. The 
fraction of the physics background (B~ — > D*°7r~) with 
respect to the signal is common to the three decays and 
the fraction of the B~ — > Dcp+n~ is common to the 
two Dcp modes. As determined from simulation, these 
modes are the only significant contributions within the 
mass range 5.17 < M < 5.60 GeV/c 2 chosen for our fit. 

Kinematic information is given by three loosely cor- 
related observables: (a) the mass M^, calculated by 
assigning the pion mass to the track from the B decay; 
(b) the momentum imbalance a, defined as 



a = 1 -ptr/PD > 

a = -(1 -PD/Ptr) < 



if Ptr<PD; 

if p tr >PD\ 



where pt r is the momentum of the track from the B 
candidate; and (c) the scalar sum of the D momentum 
and the momentum of the track from the B candidate 
{ptot — Ptr + Pd)- The above variables uniquely iden- 
tify the invariant mass Mdk evaluated with a kaon mass 
assignment to the track from the B decay, through the 
(exact) relations [23[ 
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FIG. 1: Invariant mass distributions of B~ — > Dn~ candidates for each reconstructed decay mode. The pion mass is assigned 
to the charged track from the B candidate decay vertex. The projections of the common likelihood fit (see text) are overlaid 
for each mode. 



+ 2 



+ 



Ptot 



(2 



, PtotO- 

1 mi + ' 




if a > 0; 




Using these variables, we can write Cj n = 

Pj(M D7T \a,ptot)Pj(a,ptot) and C^ m = 
Pj(dE /dx\a,ptot), where Pj is the probability den- 
sity function for decay mode j. Distributions of the 
kinematic variables for the signals are obtained from 
samples of events from the full CDF simulation, while 
for the combinatorial background they arc obtained 
from the mass sidebands of data. The shape of the mass 
distribution assigned to each signal process (B~ — > Dir~ 
and B~ — > DK~ decays) has been modeled in detail 
from a dedicated study including the effect of final 
state QED radiation [22[. The simulation results were 
tested on high-statistics data samples of D° decays, in 
order to ensure the reliability of the extraction of the 
DK~ component in the vicinity of the larger Dtt~ peak. 
Exponential functions were used to model the mass 
distribution of combinatorial background for each mode. 
The normalization and the slope of these functions are 
independently determined in the maximum likelihood 
fit. The particle identification (PID) model of the 
combinatorial background allows for pion and kaon 
components, which are free to vary in the fit. 

A large sample of D* + —> D°(—> K^tt + )tt + decays was 
used to calibrate the dE / dx response of the detector to 



kaons and pions, using the charge of the pion in the D* + 
decay to determine the identity of the D° decay products. 
The calibration includes the dependence of the shape and 
the average of the response curve on particle momentum, 
and the shape of the distribution of common-mode fluc- 
tuations. The calibrated dE / dx information provides a 
1.5 a separation power between pion and kaon particles 
of pt > 2 GeV/c. Uncertainties on the calibration pa- 
rameters are included in the final systematic uncertainty 
of A C p+ and R C p+ M- 

The B~ —¥ DK~ and B~ Dir~ signal event yields 
obtained from the fit to the data are reported in Table |U 
The fraction of the B~ — > tt + tt~K~ was set by the fit to 
its lower bound at zero, compatible with the expectation 
of a negligible contribution, and will be ignored in the 
following. The uncorrected values of the double ratio 
of branching fractions Rcp+ an( l of the CP asymmetry 
Acp+ obtained from the fit are Rcp+ = 1-27 ± 0.24 and 
A C p+ = 0.39 ± 0.17. In the fit, R C p+ and A C p+ are 
functions of the fractions [fj in Eq. [7] and the total 
number of events in each subsample. 

As a check of the goodness of the fit, and to visualize 
better the separation between signal and background, we 
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TABLE LB — > DK and B — > D-k event yields obtained from the fit to the data. 



D mode B+ -> Dn+ 


B~ -t Dtt- 


B+ -> M+ 


B" -> DJf" B+ - 


-> [h~h+]K+ B~ 


-> [/i-/i+]Jf- 


K~n + 3769 ± 68 


3763 ± 68 


250 ± 26 


266 ± 27 






K + K~ 38f±25 


399 ± 26 


22 ±8 


49 ± 11 


3±1 


3±1 


k+k~ 101 ±13 


117 ±14 


6±6 


14 ±6 







plot distributions of the relative signal likelihoods: 



RL 



pdf(B -> DK) 



pdf(B — > DK) + pdf (background) 



(8) 



where pdf(B — > _Di^) is the probability density under the 
signal hypothesis, and pdf (background) is the probability 
density under the background hypothesis (including both 
physics and combinatorial backgrounds, with their mea- 
sured relative fractions). These distributions are com- 
pared to the prediction of our fit in Fig. [2j showing a 
very good agreement. In addition, we plot projections of 
the fit on the invariant mass distributions, both for the 
entire sample (Fig. [1]), and for a kaon-enriched subsam- 
ple, where the interesting B~ — > DK~ components have 
been enhanced with respect to the B~ — > Dtt~ by means 
of a dE/dx cut (Fig. [3]). All these projections show very 
good agreement between our fit and the data. 

Some corrections are needed to convert our fit re- 
sults into measurements of the parameters of interest. 
First, we correct for small biases in the fit procedure it- 
self, as measured by repeated fits on simulated samples: 
S(R C p+) = -0.027±0.005 and 5(A C p+) = 0.015±0.003. 
These biases are independent of the true values of Acp+ 
and Rcp+ used in the simulated samples. Rcp+ does 
not need any further corrections because detector ef- 
fects cancel in the double ratio of branching fractions. 
The direct CP asymmetry Acp+ needs to be corrected 
for the different probability for K + and K~ mesons to 
interact with the tracker material. This effect is re- 
produced well by CDF II detector simulation (traced 

by GEANT 0]), which yields an estimate 

1.0178 ± 0.0023(stat) ± 0.0045(syst) E 
verified by measurements on data [2 
The corrected results are 



e(K- 

which has been 



R C p+ = 1.30±0.24(stat), 
A C p+ = 0.39 ± 0.17(stat), 



(9) 
(10) 



where Acp+ was corrected using the following equation: 



A 



N(B~ -> D° CP+ K 



CP+ 



JV(B H 



D° CP+ K+) 



N(B- 



D° CP+ K- 



+ N(B+ 



D° CP+ K+) 



(11) 



Systematic uncertainties are listed in Table [TTJ They 
were determined by generating simulated samples of 



pseudocxperiments with different underlying assump- 
tions, and checking the effect of such changes on the 
results of our measurement procedure. The dominant 
contributions are uncertainty on the dE / dx calibration 
and parametrization, uncertainty on the kinematics of 
the combinatorial background, and uncertainty on the 
physics background (B~ 



D*°tt ) mass distribution. 



Variations in the model of the combinatorial background 
included different functional forms of the mass distribu- 
tion, and alternative (a, ptot) distributions, constrained 
by comparison with real data in the mass sidebands. 

Smaller contributions are assigned for trigger efficien- 
cies, assumed B~ mass input in the fit [27j and kinematic 
properties of signal and physics background. 



TABLE II: Summary of systematic uncertainties. 



Source 



R 



cp+ Acp+ 



dE/dx model 0.056 0.030 

D*°ty mass model 0.025 0.006 

Input B~ mass to the fit 0.004 0.002 

Combinatorial background mass model 0.020 0.001 

Combinatorial background kinematics 0.100 0.020 

D-k kinematics 0.002 0.001 

DK kinematics 0.002 0.004 

D*°tt kinematics 0.004 0.002 

Fit bias 0.005 0.003 



Total (sum in quadrature) 



0.12 0.04 



In summary, we have measured the double ratio of 
CP-even to flavor eigenstate branching fractions [Eq. [5] 
R C p+ = 1.30 ± 0.24(stat) ± 0.12(syst) and the direct 
CP asymmetry [Eq. Q] A C p+ = 0.39 ± 0.17(stat) ± 
0.04(syst). These results can be combined with other 
B~ — > DK~ decay parameters to improve the determi- 
nation of the CKM angle 7. These measurements are 
performed here for the first time in hadron collisions, are 
in agreement with previous measurements from BaBar 
(R C p+ = 1.06±0.10±0.05,^ CP+ = 0.27±0.09±0.04in 
348 fb _1 of integrated luminosity @) and Belle (Rcp+ = 
1.13 ±0.16 ±0.08, A CP+ = 0.06±0.14±0.05 in 250 fb" 1 
of integrated luminosity [To| ) and have comparable un- 
certainties. 
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FIG. 2: Relative likelihood for B~ —¥ DK~ candidates for each reconstructed decay mode. The points with the error bars 
show the distribution obtained on the fitted data sample while the histograms show the distributions obtained by generating 
signal and background events directly from the total PDF of the fit composition. 




B" -» D cp+ Jt" -> [K" K + ]te" 
Negative charge: 
D™, it" 



— B 

■ B" -> D'° ji' 
B > K K* K 

combinatorial background 



4=i 



5.4 



5.5 



> 

Hi 



40 



a. 

£ 30 
c 
a 
> 

LU 

20 



1I + ^[K" K + ]3T + 

Positive charge: 
B*- 

]B^D CPt K* 
B* -> D*° Jt* 
B* -> K K* K* 
combinatorial background 



KKt: mass [GeV/c 2 ] 




KKji mass [GeV/c 2 ] 




FIG. 3: Invariant mass distributions of B~ — > Dn~ candidates for each reconstructed decay mode. The pion mass is assigned 
to the prompt track from the B decay. A requirement on the PID variable was applied to suppress the Dn component and 
favor the DK component. The projections of the likelihood fit for each mode are overlaid. The p-value for agreement of data 
with the fit is 0.95 
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